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Abstract 

The conductance through a mesoscopic system of interacting electrons coupled to two adjacent 
leads is conventionally derived via the Keldysh nonequilibrium Green's function technique, in 
the limit of noninteracting leads [see Y. Meir et ai, Phys. Rev. Lett. 68, 2512(1991)]. We extend 
the standard formalism to cater for a quantum dot system with Coulombic interactions between 
the quantum dot and the leads. The general current expression is obtained by considering the 
equation of motion of the time-ordered Green's function of the system. The nonequilibrium 
effects of the interacting leads are then incorporated by determining the contour-ordered Green's 
function over the Keldysh loop and applying Langreth's theorem. The dot-lead interactions 
significantly increase the height of the Kondo peaks in density of states of the quantum dot. 
This translates into two Kondo peaks in the spin differential conductance when the magnitude 
of the spin bias equals that of the Zeeman splitting. There also exists a plateau in the charge 
differential conductance due to the combined effect of spin bias and the Zeeman splitting. The 
low-bias conductance plateau with sharp edges is also a characteristic of the Kondo effect. The 
conductance plateau disappears for the case of asymmetric dot-lead interaction. 
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1. Introduction 

With the advancement in nanotechnology, a small number of interacting electrons can be con- 
fined in a quantum dot, to mimic the behavior of impurity atoms in a metal, thus enabling the 
study of many-body correlations between electrons. One such many-body phenomenon is the 
Kondo effect, the discovery of which in quantum dot systems |llj,|2|], has generated tremendous 
theoretical 0, El Si and experimental 10, 0] interest. Early works on the Kondo effect had 
largely focused on the effect of charge bias which applied to the two leads. Later, attentions 
have shifted to the Kondo effects of the systems including ferromagnetic leads jgl [Tol [III [l2ll 
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and the leads with spin accumulation I13H . Spin accumulation is the nonequilibrium split in the 
chemical potentials of spin-up and spin-down electrons. Accordingly, a spin bias can be realized 
experimentally, for instance, by controlling the spin accumulation at the biased contacts between 
ferromagnetic and nonmagnetic leads, and can in turn, generate a pure spin current without any 
accompanying charge current lfl4l H5I1 . Recent experimental progress has opened new possi- 
bilities for the study of spin-bias-induced transport in strongly correlated systems lfl6l \vi[ [l8ll . 
However the effect of the spin bias on the Kondo effect remains largely unexplored. More im- 
portantly, the Kondo effect in the presence of Coulombic interactions between the quantum dot 
and two adjacent leads has not been studied. 

Such dot-lead interactions warrant an investigation owing to the increasing influence of leads 
with the reduction in the lead-dot separation. The nonequilibrium (bias-driven) transport through 
a interacting mesoscopic system is usually modeled by the Keldysh nonequilibrium Green's func- 
tion (NEGF) method. The standard Keldysh NEGF method for a central quantum dot system 
coupled to two biased leads, was first developed by Meir et al. [ 19], and widely applied there- 
after m nn m m 3 . 

In this paper, we derive the general analytical expression for the current in a quantum dot 
system with interacting leads. The dot-lead Coulombic interactions are then incorporated by 
determining the contour-ordered Green's function over the Keldysh loop and applying Langreth's 
theorem. Based on the current formula, and the retarded Green's function of the quantum dot 
system, we study the low- temperature spin transport and the Kondo physics induced by a spin 
bias, in the presence of dot-lead Coulombic interactions. Two sharp peaks occur in the density 
of states in the Kondo regime. This translates into Kondo peaks occur in the spin differential 
conductance when the magnitude of the spin bias becomes equal to that of the Zeeman split of 
the quantum dot energy levels. There also exists a conductance plateau in the charge differential 
conductance at low bias, due to the combined effect of the spin bias and the Zeeman splitting. 
Since the Kondo effect primarily arises from intra-dot Coulomb interactions between electrons 
of opposite spins, the Kondo peaks and hence the position of the conductance peaks and plateau 
are largely unaffected by the strength of the lead-dot Coulombic interactions. The conductance 
plateau disappears when the symmetry between the dot-lead interactions in the left and right 
junctions is broken. 

The organization of the rest of the paper is as follows. In Sec. [2] the Hamiltonian of the sys- 
tem is introduced, and the current formula and corresponding Green's functions are derived. In 
Sec. [3] we present the results of our numerical calculation of density of states and the differential 
conductance in the presence of the dot-lead Coulombic interactions. Finally, a brief summary is 
given in Sec. [4] 



2. Model and formulation 

The Hamiltonian of the quantum dot system can be written as H — H c + H cen + Hj + H^. The 
first term is the Hamiltonian of the contacts 

while the Hamiltonian of the central region is described by: 

Hcen - 2^ E a-d cr d cr + Ud^d-\d^d[, (2) 

cr 
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where cl. {c aka -) and ct a (da) are the creation (annihilation) operators of an electron with spin 
o"(=T> 1) in the lead a(- L, R) and the quantum dot, respectively. The second term in H cen is 
the Anderson term where U refers to the on-site Coulombic repulsion, while s a kcr and e a are the 
energy levels of the two leads and the quantum dot, respectively. For simplicity, we consider 
only a single pair of levels with energies = so + crAe/2. The third term represents tunneling 
coupling between the leads and the central region: 

+ H.C.), (3) 

ktraeL,R 

where t a(T refers to the tunneling coupling constant. 

In the conventional NEGF method, the interactions in the leads are disregarded so that the 
semi-infinite leads can be described by the simple noninteracting Green's function (see Meir et 
al. Qjj] ). In this paper, however, we incorporate the Coulombic interaction between the charges 
on the quantum dot and the two leads, which may be significant given the proximity of the two 



regions. This interaction term in the Hamiltonian is given by [25] 



Hid = takdl-do-C^Cako-,, (4) 

where I a k denotes the interaction strength. We first derive the general formula of the current by 
considering the equation of motion method, and the contour-ordered Green's function over the 
Keldysh loop, as will be shown explicitly below. 

The current from the lead a through barrier to the central region can be calculated from the 
time evolution of the occupation number operator of the lead a: 

J a<T = -e{N a a) = ^Re{ V t acr G* k(r (t, t)}, (5) 
n k 

where N acr — Y,k c l : ko- Cako ' an< ^ ^ e l esser Green's function is defined as 

G^ cr (r,r') = /<4 t7 (f')^(0). (6) 
According to the equation of motion satisfied by the time-ordered Green's function G f , , namely, 

-^G' ahr (t, t') = e aktT G' ak(T {t, t') + f aa G' a {t, O 

+ Yj U(-0<T{^(f)4"( f ')^"(f')4 (T (f'))), (7) 

cr" 

where we define the central region time-ordered Green's function G^(f, t') = -/(T{J cr (f)(i^(f')}) 
with T being the time-ordering operator. Note that the Coulombic interaction between the quan- 
tum dot and the leads results in the last term in Eq. (0. This additional term prevents the closure 
of the equation of motion satisfied by G' akcr except by utilizing certain approximations. We adopt 
the Hartree-Fock approximation to deal with this term and obtain a closed set of equations. Using 
the Hartree-Fock approximation, the last term can be expressed as 

-/(r{^(f)4„(f')^"(f')cL(f'))> - ncrll (-i)(T{dAt)cl klT (t')}), 
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here = (d^d^). Thus, the equation of motion of Eq. ^ can be rewritten as 
. d 



df 



' ^akcr 



(8) 



where A ak = nl a k with n — + nj,. Thus, the Green's function G^ ka _(t, f) can be derived through 
adopting Langreth's theorem, namely 



G^Jf, t') = j dhtl^it, h)g< k(r {h,t') + G< a {t, h)g a aka (t u t')\, 



(9) 



where the time-dependent Green's function of the leads for the uncoupled system with the mod- 
ified energy level s' , are 



SahrM = if air (s' aka .)exp[-is' alar (t-t')], 
Sak^'t") = +i8(+t + t')exp[-is' aklT (t-t')], 
1-1 



(10) 
(11) 



where f a As' aka ) = { exp [(e^. - /J. acr )/k B Tj + l} . Substituting the Green's function G< k(r (t, t') 
into Eq. (0, we can obtain the formula of the current 



ie C de i „ 

+f a ,Ae + A a (e))[G' a (e + A a (e)) - G a a {e + A a (e))]}, 



(12) 



where r acr = 2nY<k \tacr\ 2 S(e - £ a far) is the linewidth function and A a (e) = nl a (e). The Green's 
function G r ^{e) are the Fourier transformations of G r <f(t) with G^it) = -i6(t)({da-(t), dl(0)}) = 
«J tr (f)l^(0))) r , G<(t) = «'(<4(0)<4-(?)> = {(d a -(t)\dl-(0))) < . For the case of proportionate coupling 
to the leads, i.e., = AFro-, the formula of the current can be simplified. Defining the current 
J = xJi - (1 - x)Jr with J a = J a i + J a ±, we can obtain the simplified formula of the charge 
current, namely 

J = d £ [f L Ae + A a ( e ))-f R Ae + A a (e))] 

rLo-(e)r So -(e) 



Yurie) + T Rcr (e) 



-ImG'Ae + AJe)) 



(13) 



Thus, finally, the current can be expressed solely in terms of the retarded Green's function. 

We now proceed to solve the retarded Green's function G^(e). The standard equation of mo- 
tion technique yields the following general relation 



= ({FuF 2 }) + (([F u H]\F 2 )Y, 



where F\ and F 2 are arbitrary operators. Note that the energy e includes a infinitesimal imaginary 
part /0 + . By considering the model Hamiltonian, the equation of motion of the Green's function 
GJJ-(f) can thus be written as 



(e - s iT ){{d <T \dl)) r = I + Yj Iak{{d lT cl kiT ,c ak Adl)Y 

akcr' 



(14) 



The above equation of motion generates three new Green's functions, i.e., {{do-c* ,,c a ka>\do)Y , 
{{d a -d-da\d a )) r and {{c a ko-\dj T }} r . We then consider the respective equations of motions of these 
Green's functions: 

(e - s afar )«Wl4» r = t a A(da\dl)Y + J] Iak((dld a ,c aklT \dl)Y. (15) 



(e - s akcr ){{d ] a .,d a ->c aktT \dl.)Y 

= ^ ta'o-'dcl^^Cako-do-'WlyY + taa-Hd^da-'d^d^Y 

a'k' 

+ ^ Co-'«^ c «'^o-'^o-l^»'" + 2 lvk{{dl„d a >>dl4a'CakAdl))'\ (16) 



(e - S(r - UW^dld^dDY 



+ ^ Co-(( d] a- d <T c a'k'o-\dl)y '+ 2^ I a>k>(( d a d l d d-cl k , cr ,Ca>k'o-'\di)Y, (17) 



a'&' a'k'cr' 



(e + Saws- - s & -s a - U)((c i a , k , g _d a -d g -\dl.)Y 



or" A:" 



+0*«^4^l4» r + J] 21**><<cl. lf rd 4tel,^c a »w]4.)y, (18) 

a ,, k ,, o Jf 

(e - S a <*' & + Sg- - £ a -)«(i -4c ff 'i:'5-l4» r 

= -(ca'k'a-dl) - U(((dg-dl + 4^o0^o-4 c Q''<t'5-l4>> r 

q,' /./ ^ ^^j/ ^/(t" d (T d 

l\d~i)Y. (19) 

To obtain a closed set of equations from the above relations, we adopt a decoupling approxima- 
tion based on the following rules JISQ: <TX> = and {{YX x X 2 \dt)Y ~ {X l X 2 ){{Y\dl)) r , 



where X and Y represent the operators of leads and the quantum dot. We shall also assume 
that higher-order spin-correlations terms can be neglected. With these approximations, Eq. ([ToT l- 
Eq. ( fT9T > simplify to 

(e - s aktT ){{dl,d a - l c oklT \dl.)Y 

= t a<T [f a ASako-)({da\ d i)Y6ao-> + <<d -4^rl4>>''<W]> (20) 

(e + s a , k , & - eg- - Sa - U){{c^ kig .d <T dg\dl.)Y 

= 04«^4^l4»'' - f a >ds a ^){{d*\dl)Y\ (21) 
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(e - s a , k , t j + eg- - So-)((do-dlc a , k ,d-\dt)y 

= t a >g\((d ( rdl.d s -\dl)Y ~ fa'd-(s a >k>d-)((do-\dl.)y 



(22) 



after ignoring all of the higher-order terms. Substituting Eq. d20l >. Eq. d2"TT > and Eq. ( |22"1 i into 
Eq. ( TTTl l. we then obtain the following: 



(23) 



In the above, the self-energies 'Lq -,\o-,2o- are defined as 

|2 



-Oo- 



Z \tg'i 



.,. £ ^a'k'cr 
a k 



n'k- 



where i = 1,2, A^,_ = 1 and A®,- = f a 'g-(e a . k >d-)- 

Subsequently, by substituting Eq. (fT5T l. Eq. (l20l > and Eq. d23l into Eq. (TBI , we obtain a relation 
involving the Green's function G^(e) = {{d tT \d] T )Y: 



(e - s lT - X 3a - ^{{do-WDY = l + (U + E 5(r )«d £r 4^l4» r 



where 



■£■30" — ^ ] Iakfacr'{ s akcr')-> ^4cr — ^ ^ 



akcr' 



£ ^akcr 



1 + 



lak facr ( ^akcr) 



-S 



(e - Safer) 2 ' 



(24) 



Finally we can obtain the analytic form of the required Green's function, namely, 

g;(6) = «4 r i4» r 

c - e (T -U - £ 0t7 - - Sip. + (£/ + £ 5o -)» g - 

(e - So- - S 3o - - S4o-)(e - fio- - U - Xqo- - s io-) + (f + ^5^)^20- 



(25) 



In the absence of any interaction between the quantum dot and the two leads, two of the self- 
energy terms reduce to zero, i.e., 1,^ = £50- = 0, while 



ak 



■'Act - M)n- 



£ £ako~ 



Eq. d25l l is the Green's function of a quantum dot system with Coulombic interactions within the 
dot, and between the dot and the two leads. For simplicity, we assume the case of strong intra-dot 
Coulomb interaction (i.e., the infinite-U limit), for which the Green's function G' a (e) reduces to: 



G;(e) = {{da-\dl)Y = 



1 - tin 



£ — En- — Sjo- ~~ ^3o- _ X4cr 
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(26) 




Figure 1 : (a) and (b) are schematic diagrams of the density of states Po-(e) for a quantum dot symmetrically coupled to 
two leads with Lorentzian linewidth of 2D. The quantum dot has two spin states with energies = -2.9, = -3.1 and 
an on-site interaction U — » oo. The linewidth is chosen to be D = 500 and the temperature is T = 0.005. The spin bias 
is V, = 0.3 and the chemical potentials are /i L j = —[III = V s /2 and /i R ^ = -)Irj = V s /2. Solid, dashed and dotted curves 
correspond to interaction parameters of It = Ir = 0, //, = /« = 0.001 and It = 0.001, Ir = 0.005, respectively. 



where the new self energy "LL is defined as 

Y> _ \ * fa' (ri^a' k'd-}\ta'd-\ (27) 

Obviously, the self-energies terms "E^ and represent the effects of the interaction between 
the quantum dot and the two leads. Note that the Green's function is dependent on the occupation 
number n a which is given by the formula - -i J(de/2K)G^.(e). The lesser Green's function 
G^(e) in the formula may be obtained by adopting certain approximation, for instance, the non- 
crossing approximation] 3] and the ansatz method[28). Alternatively, one can directly calculate 
the integral J deG^(e) via the following relation: 



C de 



ImG; rL " /Ltr(6) + W ^ (6) . (28) 

Thus, from Eqs. ( |26| | and ( |28T i. one can evaluate G r a (e) and n a self-consistently. Finally, the 
converged value of is then used to calculate the charge current via Eq. (ITU . 



3. Numerical results 

Having derived the nonequilibrium transport for the general case of a central region coupled 
to interacting leads, we investigate the spin transport through a quantum dot system with a 
spin bias V s applied at the two leads, namely ^ = -[Ili = VJ2 and fi K i = -Hr-\ = V s /2. 
Our focus is to analyze the effect of the lead-dot Coulomb interaction on the spin-transport 
properties. Firstly, we start to study density of states of the quantum dot according to the re- 
lation Po-(e) = -(l/7r)ImGJ r (e). In the following numerical calculations, we assume that the 
quantum dot symmetrically couples to two leads with Lorentzian linewidth of 2D, namely 
T Lo -(e) = T Ro -(e) = y D 2 /2(e 2 + D 2 ), with y = 1 as the unit of energy and D = 500. As 
for the dot-lead interaction, we adopt a flat-band profile, i.e., I a (e) = I a 9(D - |e|). 
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Figure 2: (a)-(b) Schematic energy diagrams of the Kondo density-of-state peaks and the electrochemical potentials of 
the two leads, for the case of Ae^ > V s and Ae^ < V s . (c)-(d) The differential charge and spin conductance G c and 
G s plotted as a function of the spin bias V s for different lead-dot interaction strengths. Solid, dashed and dotted curves 
correspond to the case of It = Ir = 0, II = Ir — 0.001 and //, = 0.001, Ir = 0.005, respectively. The energies of 
the quantum dot are chosen to be £^ = -2.9, £y = -3.1, so that Ae-j = 0.2. Other parameters are the same as those of 
Figure [T] 



As shown in Figure [TJ the spin-up and spin-down density of states are plotted in the presence 
of dot-lead Coulombic interactions described by Ii and Ir. In the absence of the dot-lead interac- 
tion, i.e., Ii — Ir = [solid line], a broad main peak is observed at e ~ -2.5, which is associated 
with the renormalized level of the quantum dot. In addition, there are two sharp Kondo peaks 
at energies e ~ and e « 0.35 [see Figure [Tfa)]. The Kondo peaks for spin <x arise from the 
contribution of the self-energy l,' 2o _, due to virtual intermediate states in which the site is occu- 
pied by an electron of opposite spin <r |3]. The real part of El grows logarithmically near the 
energies = fi aS - + As^ with Aeo- = Sg- - e Ct , due to the sharp Fermi surface at low temperature. 
This logarithmic increase translates into peaks in the density of states near those energies. Ac- 
cording to the same physical explanation, we can deduce that the Kondo peaks of the spin-down 
density of state pi should occur at energies e = -0.05 and e = -0.35, as can be confirmed from 
Figure [Tib). In the presence of the dot-lead interaction, for example II - Ir - 0.001 (dashed 
lines), the position of the broad main peak is shifted. However, the positions of the Kondo peaks 
are not affected by the dot-lead interaction which induces the self-energies Y,^ lT and and do 
not contribute to the Kondo effect. When the interaction strength is increased, e.g., Ir = 0.005 
(dotted line), the position of the broad main peak is shifted further to higher energy compared to 
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Figure 3: The differential charge and spin conductance G c (a) and G s (b) versus the spin bias V s , for temperatures 
T = 0.005 (solid line), T = 0.05 (dashed line) and T = 0,1 (dotted line). The lead-dot interaction strength is set at 
II = Ir = 0.001. Other parameters are the same as those of Figure [TJ 



the previous case, but the positions of the two Kondo peaks remain invariant. 

Next, we investigate the effect of the dot-lead interaction on the charge and spin differential 
conductances, which are defined as G c - dJ/dV s and G s - dJ s /dV s , where J s = h{J-^ - J\)j2e, 
respectively. The differential conductance G c and G s are plotted as a function of the spin bias 
V s for different interaction strengths as shown in Figure To explain the observed trends in 
the spin and charge differential conductance, we sketch the electrochemical potentials in the two 
leads, and superimpose on it the Kondo peaks in the density-of-states [see Figs. [3a) and (b)]. 
We observe a plateau in G c over the bias interval of \V S \ < 0.2, for the cases of Ii — Ir — and 
Il-Ir- 0.001. However, the conductance plateau is destroyed in the case of asymmetric lead- 
dot interaction, i.e. Il^ Ir, and the charge differential conductance assumes a negative value over 
the same bias interval [see Figure [2c)]. The spin differential conductance G s shows two Kondo 
peaks at \V S \ = 0.2, irrespective of the symmetry or strength of the lead-dot interactions [shown 
in Figure |2jd)]. G s is also enhanced with increasing strength of the dot-lead interactions. We 
note that the spin-up and spin-down electrons flow along opposite directions, i.e., J = \Jj\ - |7J 
and J s oc + \J^\, and that the two currents 7| and Ji can be different due to the energy 
splitting Asf. To explain the above conductance dependence on V s , we refer to the schematic 
diagram of the Kondo density-of-state peaks [see Figure |2ja)-(b)]. The Kondo peaks begin to 
enter the spin-bias conduction window when the spin bias is increased beyond \V S \ = Ae-f. This 
results in the spin differential conductance G s having two Kondo peaks at V s = ±Asf. The 
entry of the two Kondo peaks into the conduction window also reduces the difference in the 
magnitude of /j and /|, thus resulting in a sharp drop (plateau step) in the charge conductance at 
V s = +Asf . The conductance plateau can thus be attributed to the combined effect of the spin bias 
in the leads and the Zeeman splitting in the QD. The dot-lead interaction tends to increase the 
coupling between the leads and the QD, thus resulting in a general increase in the charge and spin 
differential conductances. In the presence of asymmetrical dot-lead interaction, i.e., Ii = 0.001 
and Ir = 0.005, the symmetry in the transport across the QD is broken, and thus the conductance 
plateau disappears. Two conductance dips occur at |V,| < Aej due to the contribution from the 
Kondo peaks in the density-of-states. 

Finally we study the temperature dependence of the differential conductances. As shown 
in Figure [3] the two Kondo peaks in G s at V s = ±0.2, and the conductance plateau in G c in 
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the bias interval \V S \ < 0.2 can be clearly observed at a low temperature of T — 0.005. With 
increasing temperature, e.g., at T = 0.05, the Kondo peaks become thermally broadened, while 
the plateau in G c sharpens into a peak profile. With a further increase in temperature to T = 0.1, 
the plateau in G c is almost completely suppressed. These changes can be largely attributed to the 
thermal distribution of electrons about the electrochemical potential in the leads. The thermal 
distribution in turn affects the self-energy El of the intra-dot Coulomb interaction, which is 
primarily responsible for the Kondo resonances in the density-of-states [see Eq. d27H . 

4. Summary 

In this work, we analyze the spin-transport properties of a quantum dot system driven by spin 
bias in the presence of dot-lead Coulombic interactions. The transport property is discussed on 
the basis of Keldysh nonequilibrium Green's function framework. According to the equation- 
of-motion technique and Langreth's theorem, we derive the analytical expression of the current 
through the quantum dot in the presence of the dot-lead Coulombic interaction. Our numerical 
results show that although the interaction can renormalize the energy levels of the quantum dot, 
they leave the position of the Kondo peaks in the density of states unchanged. This is because the 
Kondo effect arises primarily for intra-dot Coulomb interactions involving electrons of opposite 
spins. The Kondo resonances in the density of states translate into peaks in the spin differential 
conductance when the magnitude of the spin bias is equal to that of the Zeeman energy split in 
the quantum dot. There also exists a plateau in the charge differential conductance at low bias, 
due to the combined effect of spin bias and the Zeeman energy splitting. The position of the 
steps of the conductance plateau can also be attributed to the Kondo effect. The strength of the 
Coulombic lead-dot interactions affects the magnitude of both the spin and charge conductances. 
Furthermore, in the presence of asymmetrical dot-lead interaction strengths, the plateau in the 
charge conductance disappears, and is replaced by conductance peaks. Finally, the temperature 
dependence of the differential conductances is qualitatively discussed. 
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